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Abstract 

A Pixelized Photon Detector (PPD) is a generic name for the semiconductor devices operated in the Geiger-mode, such 
as Silicon PhotoMultiplier and Multi-Pixel Photon Counter, which has high photon counting capability. While the 
internal mechanisms of the PPD have been intensively studied in recent years, the existing models do not include the 
avalanche process. We have simulated the multiplication and quenching of the avalanche process and have succeeded 
in reproducing the output waveform of the PPD. Furthermore our model predicts the existence of dead-time in the 
PPD which has never been numerically predicted. For serching the dead-time, we also have developed waveform 
analysis method using deconvolution which has the potential to distinguish neibouring pulses precisely. In this paper, 
we discuss our improved model and waveform analysis method. 
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1. Introduction 

A Pixelized Photon Detector (PPD) |l| is a novel 
semiconductor photon sensor which has a single pho- 
ton level sensitivity and a linear relationship be- 
tween the gain and the operation voltage [2| . More- 
over the output waveform of the PPD has two com- 
ponent and is not described by the typical single 
exponent function 0, 0] ■ While the internal mecha- 
nisms of the PPD have been intensively studied in 
recent years, the existing models incorporate arti- 
ficial manipulation and are nothing more than the 
phenomenology which explains above characteris- 
tics. Thus, we have developed a new model based on 
semiconductor physics and succeed in reproduction 
of the all characteristics. 
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Fig. 1. A simple circuit of the PPD. Each pixel consists of a 
diode and a quenching resistor. 

A traditional model is shown in Figure [T] where 
V op is the operation voltage, Cd the capacitance 
value of the diode and R q the quenching resistor 
value. The main difference between our model and 
the previous ones is the precondition of the applied 
voltage to the diode-part in the PPD (Vd). The ex- 
isting models fixed the minimum value of Vd at the 
breakdown voltage (Vb r ) and the avalanche multi- 
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plication is forced to be terminated when Vd reaches 
Vbr- On the contrary, we did not request such a con- 
straint and the avalanche process is terminated due 
to the internal physics: charge transportation, im- 
pact ionization and circuit equation [H, 0| . On the 
other hand, Vd may drops below Vbr, ancl hence the 
gain linearity is not trivially determined. 

We have though calculated and we have in fact 
confirmed the gain linearity and the difference be- 
tween V m in and Vbr, therefore the existence of sig- 
nificant dead time of the PPD is predicted. 

In this paper, we firstly introduce our model of 
the PPD and the prediction of the dead time. Next, 
we report the waveform analysis method which is 
effective to search the dead time. 

2. A new model incorporating the avalanche 
process 

Figure shows our model and characterized by 
the charge of created carriers q(t) in the diode which 
is time dependent in order to treat the transient 
multiplication correctly [3]. In this Figure, I(t) 
means the observable current and C q describes a 
capacitance located between the diode-part and the 
resistor-part, which has already reported by ITC- 
irst group 

UMM- Note that, V d (t) is determined by 

both q(t) and I(t): V d (t) = V op - dt _ Us _ 

ing this formula, the circuit equation is calculated 
as follows: 



dV d (t) 
dt 



R q {Cd + C q ) 



(Vop - V d {t) 



^^ r )-(i) 



The variation of the created carriers ^4P- is de- 

dt 

scribed by the densities of electrons and holes in the 
diode (p e , ph ) , the ionization rates (a e , a;, ) , the drift 
velocities (v e , Vh) and is expressed by: 



dg(t) 
dt 



E 



Pi{x)^(v d{t) )Vt dx 3 , 



(2) 



The ionization rate strongly depends on the electric 
field (E(x)) and is proportional to E(x) exp(— -g^y), 
where B is a constant. We caluculated the E(x) of 
the typical PPD [n| and obtained the almost uni- 
form E{x) in the multiplication layer in the diode. 
Thus, we assumed that a e and ah are the fuction of 
V d . 

Equations (TT]), ([2]) describe the evolution in the 
PPD and can be sequentially solved. Consequencely, 
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Fig. 2. A new equivalent circuit proposed in this paper. 

the observable current is culculated from the follow- 
ing equation: 



/(*) - c d 



dV d (t) , dq(t) 



dt 



dt 



(3) 



Figure [3] shows that the waveforms expected from 
equation ([3]) (solid) agree well with the measured 
ones (dashed). Note that, the transmission proper- 
ties of the measurement system are already consid- 
ered. 
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Fig. 3. The waveform expected from our proposed model 
and the measured waveform. 
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Figure 2] shows the time variation of Vd and indi- 
cate that Vd drops below Vb r , therefore a few nano 
seconds dead-time is expected. We ignored the the 
carrier lifetimes and recombinations which decrease 
the number of the carriers. If the effects are taken 
into account, the dead-time seems to be less than 
predicted, however it will still be significant. 
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Fig. 4. The waveform expected from our proposed model 
and the measured waveform. 



3. Waveform analysis 

For the precise measurement of average waveform 
of PPD, or for the study of crosstalk, afterpulsing, 
recovery, etc., it is effective to use and analyze wave- 
forms. We are using a digital oscilloscope and mea- 
suring the waveform for studying basic characteris- 
tics of PPDs. Consequently, waveform analysis al- 
gorithm is important. Here we introduce the decon- 
volution method we are developing. 

If the response function h(t) to a particular im- 
pulse is known, deconvolution using Fourier analysis 
is effective to detect pulses precisely for linear re- 
sponse systems. Here we describe the outline of this 
algorithm. 

Suppose that the observed waveform is written as 
follows: 

K«)-/*T*e-«««-o + «e) 



= / dt'd(t')h(t - t') + n(t) 



(4) 



where d(t) = J2j djS(t — tj) denotes the array of 
impulses, which is referred as comb function, and 
n(t) denotes the noise component of the observed 
waveform. Following this assumption, out problem is 
to estimate the comb function as precise as possible 
and to find and measure the position and the size of 



peaks {a.,-, tj} from the estimated comb function. In 
the frequency domain, Eqf4]is represented as 



Y(w) = D(w)H(w) + N(lu) , 



(5) 



where Y(lo), D(u>), and H{uS) are Fourier transfor- 
mations of y(t), d(t), and h(t), respectively. Consid- 
ering the ideal case that the noise component is ne- 
glected, Eq|5]is replaced to 



Y{lo) = D(uj)H(uj) 



(0) 



Dividing Y(uj) by H(ui), D(lo) is obtained, and the 
problem is solved as this is equivalent to the comb 
function d(t): 



d{t)=T 



1 



H{ui) 
(h D ®y){t) 



Y(co) 



= T- 1 [H d (u)Y(uj)] 



(7) 



where Hd(ui) = if(u;) _1 and the symbol (£> denotes 
the convolution. 

Next we cope with the real calse. In the real case 
that the noise component exists, the deconvolution 
method cannot be used directly because the power 
spectrum of the response function |7J(o;)| 2 is typi- 
cally large in lower frequency region and is small in 
higher frequency region. Consequently the deconvo- 
lution filter enhances the high frequency region of 
the observed waveform y(t). As the power spectrum 
of the noise component is roughly flat in frequency, 
this manner of the deconvolution filter emphasizes 
the noise component, and as a result the deconvo- 
lved waveform oscillates keenly. This is not the de- 
sired one. 

To avoid this problem, couple of the deconvolu- 
tion filter with some low-pass filters (LPFs) is effec- 
tive. Several types of LPFs arc selectable e.g. Wiener 
optimal filter, at present we adopt the proper LPF: 



LPF{lo) 



I 



y/1 + (W/W C ) S 



(8) 



where uj c is the cutoff frequency that is determined 
from the relation of the time resolution and the con- 
fidence level of the peak detection. The higher the 
cutoff frequency the lower the SNR, and the lower 
the lower the time resolution. Using this deconvolu- 
tion method we can split pile- up neighboring pulses 
up to 2 [nsec] intervals. Figure [5] shows the scheme 
of the deconvolution method. 
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Raw waveform 
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Extracting Peak Information 



Peak ID Time [nsl Amplitude Width [ns] 

200.35 8.64x10 2.35 

1 227.23 7.90 xlO 2 2.21 

2 239.70 8.21 xlO 2 2.22 

3 831.92 8.52 xlO" 2 2.24 



Fig. 5. The scheme of the waveform analysis using deconvo- 
lution method. 

4. Conclusion and future prospects 

We have developed a new model of Pixelized Pho- 
ton detectors and have succeeded in reproducting all 
characteristics, such as spike component in the wave- 
forms and linearity relationship between the gain 
and the bias voltage. Furthermore we our model pre- 
dicts the existence of dead-time in the PPD which 
has not been numerically predicted. 

Our waveform analysis method has the poten- 
tial to distinguish neighbouring pulses upto about 
2[nsec] interval and to determine the amplitude of 
each pulse precisely. Applying this method, various 
precise analyses are possible, e.g. measuring the re- 
covery process to search dead-time, or inspecting the 
frequency distribution of afterpulsing, etc. 

Now, the measurement of the dead-time and the 
verification of our model of PPD is on going. We are 
also challenging to study the characteristics of PPDs 
using TCAD semiconductor simulation technology. 
This study will lead to strategic developments of 



PPDs, especially in developing optimized PPDs for 
applied experiments in the future. 
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